


MATHEMATICS MAGAZINE 


Formerly NATIONAL MATHEMATICS MaGazinz, founded by S. T. Sanders 


EDITORIAL STAFF 


EXeEcuTIvE ComMit1TEE: D. H. Hyers, A. D. Michal, Glenn James. 
Current Papers AND Books: H. V. Craig. 

HisToRY AND HUMANISM: G. Waldo Dunnington, A. W. Richeson. 
TEACHING OF MATHEMATICS: Joseph Seidlin, L. J. Adams, C. N. Shuster. 
PROBLEMS AND Questions: C. J. Jaeger, H. J. Hamilten. 
MATHEMATICAL MISCELLANY: Marion E. Stark. 


GENERAL: W. E. Byrne, Irby C. Nichols, C. D. Smith, H. L. Smith, P. K. Smith, 
S. T. Sanders (Emeritus). 


Address general editorial correspondence to the managing editor, GLENN JAMEs, 
University of California, Los Angeles 24, California, and special papers to the editors 
of the departments for which they are intended. 


Manuscripts should be typed on 8$’’x11” paper, double-spaced with 1’’ margins. 
Because of the wide variety of papers published in this magazine, it is highly desirable 
that each introduction should state what the paper accomplishes and that each paper 
should be accompanied by a general statement of what mathematical training is re- 
quired to read it. Expository papers which can be read by men who are not special- 
izing in the field in which a given paper is written, are especially desired. 


The MatHEematics MaGazineE is published by the Managing Editor, bi-monthly 
except July-August. Ordinary subscriptions are $3.00, sponsoring subscriptions $10.00 
single copies 65 cents. A Subscription card enelosed in your magazine is an invitation 
to subscribe, to renew your subscription or to pay past pledges, as the case may be. 
Send subscriptions to GLENN JAMES, 14068 Van Nuys Blvd., Pacoima, California. 
Application for entry as second-class matter is pending. 


SPONSORING SUBSCRIBERS 


E. F. BECKENBACH RosBert B. HERRER S. T. SANDERS 

H. W. BECKER ALFRED HuME Cart N. SHUSTER 
Currrorp BELL Donacp H. HyYers H. A. Simmons 
Homer V. Craic I-tinots State Norma, 4 I. S. SOKOLNIKOFF 
Josern W. CRerLey GLENN JAMES D. Victor STEED 
Henry F. Ferris CLypbE E. Love E. M. TINGLEY 
Harry M. GEHMAN ARISTOTLE D. MICHAL H. S. VANDIVER 

J. GinssurG BAYARD ROBINSON ALAN WAYNE 
Tuomas GRISELLE Ira ROSENBAUM W. L. G. WitiaMs 
Gorpon HALL N. E. Rutr 


Concerning everything pertaining to advertising in this magazine, contact HomER 
V. CraiG, University of Texas, Austin 12, Texas. 





CONTENTS 


The Operational Calculus—by F. D. MURNAGHAN- 


Anyone who has had a course in Advanced Calculus including elementary 
differential equation theory and elementary complex function theory should 
be able to read this paper on Operational Calculus. The reader must have some 
acquaintance with matrices but no deep knowledge of matrix theory is neces- 


sary. 
A Novel Algorithm at the Freshman Level—by GERALD B. HuFF 
An understanding of undergraduate mathematics through the calculus 


should be enough for reading this paper Possibly it touches on a point or so that 
might be unfamiliar to such a student, but that is the level of the ideas. 


Current Papers and Books 


History aud Humanism 


A French Mathematician of the Sixteenth century—by V. T. THEBAULT. 
Teaching of Mathematics 

Proofs of the Fundamental Theorems of Spherical Trigonometry—by C. 
CLARK. 
Problems and Questions. 


Mathematical Miscellany -_- 


Our Contributors. : Cover 


The next issue of MATHEMATICS MAGAZINE will contain, in part, the first install- 
ment of a semi-expository paper on ‘The Generalized Weirstrass Approximation 
Theorem” by Marshall H. Stone, a paper on ‘‘Congruence Methods as Applied to 
Diaphantive Analysis” by H. S. Vandiver, and the first of a series of articles telling 
in simple language ‘‘the story’’ of various mathematical courses designed to bridge the 
chasm between the general reader and the abstract world of mathematics. 











The Operational Calculus 
by F. D. MurRNAGHAN 


It is now generally recognized that the operational calculus fur- 
nishes a powerful method for obtaining particular solutions of a system 
of linear differential equations with constant coefficients. Most expo- 
sitions of operational calculus, however, do not make it quite clear 
just what particular solutions of the system of linear differential equa- 
tions are furnished by the operational calculus and many are restricted 
to the case of a single differential equation. We think, then, that an 
expository article on the operational calculus may serve a useful purpose. 
The principal novelty in this exposition lies in the fact that we do not 
confine ourselves to the case of a normal system of linear differential 
equations, i.e., one in which the determinant of the coefficients of the 
highest order derivative which appears in the system is different from 
zero. We treat also the abnormal case in which this determinant is 
zero. In the normal case the derivatives of order lower than the highest 
which occurs in the system may be assigned arbitrary initial values 
while in the abnormal case these lower order derivatives must satisfy 
certain relations. 


1. The Laplace Transformation. Operational calculus is based on 
an operation (known as the Laplace Transformation) which sends com- 
plex valued functions of a real variable ¢ into functions of a complex 
variable s. The functions of the real variable ¢ are defined only for 
120, i.e., for non-negative values of ¢, but it is convenient to extend 
their range of definition by the convention that they are zero when t<0. 
It is sufficient for our purpose to confine our attention to piecewise- 
continuous functions. If f(t) isone of these piecewise-continuous functions 
we denote by f(+0) the limit at ‘=0 of f(t), it being understood that ¢ 
tends to zero through positive values; thus if f(¢) is continuous at ¢=0, 
f( +0) =f(0) =0, but if f() is not continuous at ¢=0, f(+0) is not neces- 
sarily equal to f(0). We suppose that there exists a complex number 
Zy=Xo+ ot such that the improper integral 

@ 
e~ *0'f(t)dt 
0 
exists. Then if s=x+ 7 is any complex number which lies to the right 
of go (i.e., for which x >x») the improper integral 


co 


¢(s) = e~?'f(t)dt 
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exists and is analytic at z.* We term this integral the Laplace Trans- 
form of f(t) and we denote it by L(f).. Thus 


If L(f) as defined at s=%, L(f) is analytic over the half-plane 
R(s) > R( so) where R(s) denotes the real part of s. 


It follows from the very definition that L(f) is a linear operator, i.e. 
Licf) =cL(f);¢ any constant; 
Lif, + fo) = L(fi) + Life & 


If we denote by 1 the unit-function, i.e., the function which is 1 when, 
i>0 and 0 when ¢<0 (the value at /=0 being unimportant) we have 


| 
\ L(l)=—; R(z)>0 


~ 


and if ¢ is any complex number we have 





| 
L(e"*) = - R(z+c)>0. 
é re ) 


It follows readily that 


z ’ a 
LOO Gl) tem L(sin at) =——— ;  R(s)>|I(a 


9 
- a~ 


2*-+a? 27+ 
where J(a) denotes the imaginary part of the (complex) number a. 
Furthermore if p is any positive integer 
L | tPf(t)} =(—1)6°L(f) 
where 6 denotes differentiation with respect to s. For example 
2az 
(2?-+-a?)? 


Of particular importance is the following result: 


L(t sin at) = ; R(s)>|I(a)]. 


If L(f) is a proper rational fraction this proper rational fraction has 





near s= @, the power series development 
. f(+0) . Df(+0) D"~'f(+0) 
Lp=————+ ie neo... 


where D denotes differentiation with respect to t and Df(+0), for example, 
is the limit at t=0 of Df, it being understood that t-+0 through positive 
values. 


*The proofs of the statements made in Sections 1, 2 and 5 of this article may be 
found in the chapter entitled Operational Calculus of the book Applied Mathematics 
by the author, which is now in press, 
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On the other hand a rational fraction which is not proper is not 
the Laplace Transform of any function of ¢ (since |Z(f)| may be made 
arbitrarily small by merely making R(z) sufficiently large). 
Finally we state the following uniqueness theorem: 


f(t) is determined unambiguously by its Laplace transform L(f) at any 
point t at which f(t) 1s continuous. 
This fundamental result is a consequence of the fact that L(f) is, 


when regarded as a function of y, the Fourier transform of the function 
of t which is zero when t<0 and which = (27)te~’‘f(t) when ¢>0.. 


2. 7h coivo'ution process. 
If f and g are any two piecewise-continuous functions of / the integral 
t 
f(r) g(t—r)dr 
0 
is a continuous function of t which we denote by f * g. It is easy to show 
that f *,=g * f and we term f * g the convolution of f and g._ If Df exists 
and is piecewise-continuous then, at any point ¢ where g is continuous, 
/* ¢ is differentiable, its derivative being furnished by the formula 


D(f * g)=f(+0)g+ (Df * g). 
The convolution process is associative, 1.¢., 
(f * 9) *h=f * (g* h) 


and so we denote the common value of these two expressions by the 
svmbol f *g*h. The significance of the convolution process for opera- 
tional calculus lies in the following theorem: 


If L(f) and L(g) are both defined at 2 the Laplace Transform of f * g 
exists and has the value L(f) L(g) at any point z for which R(z—2%9) >0. 





t 
Since f*1= f(r)dr and since L( ws , R(z)>0, 
‘ 2 
t . 
we have Lyf findre= scl 
\ Jo 


it being understood that the real part of s is positive and that z lies to the 
right of a point 2 at which L(f) is defined. This result may be conven- 


iently phrased as follows: 
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Integration over [0,¢] is reflected in the domain of Laplace Transforms 
by division by z. 


It readily follows that if f(t) possesses a piecewise-continuous de- 
rivative D(f) and if D(f) possesses a Laplace Transform at z= 2» then at 
any point z for which R(s—2)) >0, L(f) exists and 


L(Df) =2} L(f)} —f(+0). 
In particular if f(+0)=0 we have L( Df) =sL(f); in other words 


Differentiation with respect to t is reflected (for those functions for 
which f(+0)=0) in the domain of Laplace Transforms by multiplication 
bys. If f(+0)£0 this multiplication by z must be followed by the adjust- 
ment involved in subtracting f(+0) from z L(f). 


We conclude with the following extension of this result: 

Let f=f(t) possess a piecewise-continuous mth order derivative 
which possesses a Laplace Transform at s=%); then L(f) exists at any 
point z for which R(z—z») >0 and 


L(D"f) = 3"L(f) —2"-f( +0) —2""*Df(+0) —--- —D"- (+0). 


The particular case n=2 is of especial importance in applications of the 
theory; for this case we have 


L(D*f) =2*L(f) —z f(+0) —Df(+0). 


3. The operational solution of a single linear differential equation 
with constant coefficients. \We write the differential equation, of order 1, 
as follows: 


P(D)y=f. 


Here D denotes differentiation with respect to the independent variable / 
and P(D) is the polynomial operator. 


P(D) =coD*"+e,D""!+- - ++ +n. 


where ¢o,:++,€, are constants. f=f(t) is a given piecewise-continuous 
function of ¢ which is zero if t<0 and which possesses at some point 
Zo=Xotyot a Laplace Transform. The problem is to determine y= (/) 
so as to satisfy the differential equation and to possess, together with its 
derivatives up to the (n—1)st, inclusive, prescribed initial values »(+ 9), 
Dy(+0),---, D™-y(+0). To solve this problem we assume that y has 
a piecewise-continuous mth order derivative D”y which possesses at some 
point z, whose real part is greater than xp, a Laplace Transform. We 
have, then, the series of relations 
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L(Dy) =< L(y) — (+0) 
L(D*y) = 2*L(y) —zy(+0) —Dy(+0) 


L(D"y) = 2"L(y) —2"~1y9(+0) —--- -—D""!y(+0) 
and these imply that 
L{ P(D)y} =P(2)L(y) — (Coa + +» -+6n-1)9(+0) — «+» —coD"~1y(+0). 


Choosing s so that it lies to the right of every zero of the polynomial 
P(s) we are sure that P(z)#0 and so 


; L(f) co" t+ rade “+Cn 1 
Liy)s—— 
’ P(z) o P(z) 





Co 
a D"1y(+0) 
(3) 


(since P(D)y=f). The initial conditions expressed by the statement 
that y and all its derivatives up to the (7—1)st, inclusive, are zero at 
t= +0 yield a particularly simple formula for L(y), namely, 
1) x2. 
P(z) 
We shall term the solution of the problem which corresponds to these 
particularly simple initial conditions the rest solution and we shall de- 


note the cest solution by r=r(t). Thus 


Lif) 


L(r) = 
7 P(z) 


If so(t), s,(2),-+ +, S,-1(¢) are the functions of ¢ whose Laplace Transforms 


cos" 1+ le de +Cn-1 Co 
are ’ = ’ 
P(z) P(z) 





respectively, we have 

L(y) = L(r)+9(+0)L(s0) +--+ -+D"-1y(+0)L(s,_1). 
Since any function of ¢ which possesses a Laplace Transform is unam- 
biguously determined, at any point where it is continuous, by this 
Laplace Transform, it follows that 


y=r+y(+0)50+ - >> +D"19(+ 0) 5n-1. 


Thus y=y(t) is determined as soon as we have found the rest solution 


r=r(t) and the auxiliary functions (s9,---, Sn-1). To facilitate the de- 
termination of r and of these auxiliary functions we consider the problem 
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of determining the solution g=q(t) of the associated homogeneous ditler- 
ential equation 


P(D)u=0 
which is zero together with its derivatives up to the (7—2)nd, inclusive, 


at t= +0, while 


| 
Sa | +() =—. 
an 


The same argument that furnished L(y) yields 


P(s)L(q) =coD"~1q(+0) 


— 


I 
P(s) 





so that L(q)= 


Lif 
Since Lir) = a. =L(q)L(f) we have 


rag * f 


(for the convolution of two functions is reflected in the domain of Laplace 
Transforms by simple multiplication). Furthermore since g(+0)=0 
we have 





L(Dq) =s L(q)= 








P(s) 
and, similarly, 
2* rs 
L(D*q) = Po »- °°, L(D*- 19) = PG) 
Hence So=CoD"— "g++ +> +en-19 3 


$1 =CoD"~ "Qt + ++ +En-29 | 


Sn-1 = Coq. 


Thus the solution of the problem is furnished, in terms of the function 
g=q(t), by the formula 


y=rt+(coD" G+ + >> +6n-1)( +0) + - - > +e0gD"- *y(+ 0) 


where r=q * f. 
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The function g=q(t) is easily obtained by analysing 


| 
P(s) 





into simple fractions. In fact the function which has 


I 


- ae 


as its Laplace Transform is e*'; the function which has 


I 
(s—a)? 
as its Laplace Transform is 
p-1 
Sa en! 
(p-1)! 
] 


the function which has 





(s—a)?+b? 


] , 
as its Laplace Transform is “s e*' sin bt 
) 


° 
~~ 


(s—a)?+b? 





and the function which has 


as its Laplace Transform is {cos bt+(a/b)sin bt}e and so on. 


In parti- 
cular if all the zeros (a,-+-, a@,) of P(z) are simple we have 





P(z) =co(3—a1)° - -(3—an) 
] = ] 1 
and ee 
P(z) uu bP (ax) (Z—ax) 
n erk! 


so that git)= >, 


k 1 6P(a,) 
We shall denote by R(t) the cest solution obtained when f(t) is the 


unit function 1. Since 


oo 1 
Lil)=—. 


~ 
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] 
‘e have L(R) =———_ 
we have =P(s) 


t 


Hence R(t) =q* 1= q(r)dr. 


0 


However the expression for R(t) may best be found directly without 
first determining g(t). All we have to do is to analyse 
| 


2P(z) 





into simple fractions. In the particular case where P(s) has n simple 
zerOs (a1,° °°, @,), none of which 1s sero, we have 








l l . l 
7 4 _ 
2P(z) P(0)z » a,6P(a,)(S—ax) 
l - gat 


and so R(t) 


—————-+}- ——_—_—__—_—<—., 
P(O) 2d ap 6P (cy) 


This result is known as the Heaviside Expansion Formula; pay particu- 
lar attention to the fact that the solution it furnishes is the rest solution, 
i.e., the solution which is zero together with all its derivatives up to 
the (n—1)st, inclusive, at /=+0. 


i L(pyo ners I =<: L(f) L(R 
Since Ar) = Pis) =sL(f) —— (f) LUN) 





and since z L(f) =L(Df)+/(+0) we have 
L(r)=L(Df) L(R)+f(+0) LCR). 
Hence, since convolution is reflected in the domain of Laplace Trans- 
forms by simple multiplication, 
r=f(+0)R+(Df * R) 
t 


1.€., r(t) =f{+0)R(t)+ Df(r)R(t—71)dr. 


0 


This important result is known as the Maxwell-Boltzmann-Hopkinson 
Principle of Superposition. It furnishes a convenient method of de- 
termining the rest solution r(t) corresponding to a given applied impulse 
f(t) when the rest solution R(t) corresponding to the unit applied im- 
pulse is known, 
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Example. (D?+n*)y=A cos wt. 
Here P(z) =2?+n?. if nZ£0, q(t) =(1/n)sin nt while, if n=0, g(t) =¢. 
Case l. nO. 
Since co=1, c,=0, co=n? we have 
So=CoDg+ci=cos nt; s;=cog=(1/n)sin nt. 


The determination of r(¢) depends on whether w?2= 7? or not. 
| 


Case 1*.  w?n?, 














As A ( Z 3 
ee ee rt et eraeaneen 
(s?+-n?) (s?+w") n*=—w? s*+u?" s°+n? 
A ( ' 
so that r=- — cos wl—cos nt}. 
n*>—" 
Case 1°. w?=n? 
Az a. ( 1 ) 
Ss bearer eae tieadlias 3 apr ape 
(s?-+ n°)? 2 s°-+n? |} 
A . 
Hence r=—— tsin nt. 
2n 


Thus the solution in Case 1 is as follows: 


A 


Case 1", y(t) = ———(cos wt—cos nt)+y¥(+0)cos nt 


n*—w- 


— sin nt: 


Dy(+0) 
——— 














n 
; b A ; Dy(+0) 
Case 1”. y(t) —— tsin nt+ y(+0)cos nt-+————- sin nt. 
2 n 
Case 2. n=0. 
Since q(t) =t we have sp=1, 5, =4. 
Case 2°, wH0 
I A | Z 
L(r)= ss ng aot —_ -—- —— 
2(3*+w’") w* g 3-+w* 
so that r= — (1—cos wf). 
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‘ ei 
Case 2”. w=0 


A t? 








L(r)= so that r= 


~ _ 


Thus the solution in Case 2 is as follows: 


(1—cos wt) +9(+0)+Dv(+0)t ; 





Case 2, y(t) =—— 


w? 


At? 
Case 2”. y(t) ee +0)+Dy(+0)t. 


4. The operational solution of a system of linear differential equations 
with constant coefficients. We shall regard a system of m differential 
equations for m unknowns as a vector differentiai equation (the m un- 
knowns y!, y?,---, y" being the coordinates of the unknown vector 
y=y(t)). The vector differential equation may be written in the form 


P(D)y=f 


where, now P(D) is an mxm matrix whose elements P,*(D) (P,(D)) 
being the element in the sth row and rth column) are given polynomial 
differential operators. If the highest order derivative that appears 
in any of the m? elements P,*(D) of the matrix P(D) is the nth we may 
write 


P(D) = CyD* + C,D*"'+ ---+C, 


where the coefficients (Co,Cj---,C,) are, now, constant mxm matrices. 
Co is not the zero matrix but it may very well be singular, i.e., its de- 
terminant may be zero. We suppose that the given vector function 
f(t) possesses at some point % a Laplace Transform and that our vector 
differential equation possesses a solution y=y(¢) having a piecewise- 
continuous th order derivative which possesses, if the real part of s is 
sufficiently large, a Laplace Transform. (We understand by the Laplace 
Transform of a vector, or of a matrix, the vector, or matrix, whose 
coordinates are the Laplace Transforms of the corresponding coordi- 
nates of the given vector, or matrix.) Proceeding in exactly the sam 
way as in the case of a single differential equation, we obtain the series 


of relations: 
L(Dy) =2L(y) —y(+0) 
L(D*y) =2*L(y) —3 »(+0) —Dy(+0) ; 


L(D"y) =3"(Ly) —2""1y(+0) —--- —D"-1y(+0) 
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and these imply that 
L\ P(D)y} = P(s) L(y) — (Cos E+ + Cy) (+0) — - -» — CoD" -y(+0). 
This relation reduces to 
L} P(D)y} = P(z)L(y) 
if yi +0), Dy(+0),---, D"~y(+0) satisfy the following relations: 


Cov( +0) =0 ; 
CoD y( +0) + Cy +0) =( i 
CoD"—!y( +0) + - >> +Cr-19(+0) =0. 


If Cy is non-singular these relations yield, one after the other, the re- 
lations 


y(+0)=0; Dy(+0)=0; ---; D"-y(+0) =0. 


Whether or not Cp is singular we shall term any solution of our vector 
differential equation whose initial values satisfy the relations: 


Cov( +0) =0;--- 3 CoD"~!y(+0)+ - > +C,-19( +0) =0 


’ 


a rest solution and we shall denote a rest solution by r=r(/). Then 


P(z)L(r) =L{ P(D)r} =L(f). 


Choosing z so that it lies to the right of every zero of the polynomial 
det P(z) we are sure that P(s) is non-singular. Denoting its reciprocal 
by N(s) we have 


L(r)=N(s)L(f). 


We shall say that our vector differential equation is nermal when Cy is 
non-singular and abnormal when Coy is singular (so that det Co=0). 
When the vector differential equation is normal det P(s) is a polynomial 
in s of degree mm and each cofactor of P(s) is a polynomial in s of degree 
not greater than n(m—1). Hence each element of N(s) is a proper 
rational fraction, When the vector differential equation is abnormal 
this may not be the case. We shall consider only the case where 


l 
Each element of — N(s) is a proper rational fraction. 
In order to get some idea as to how restrictive this assumption is, con- 
sider for a moment the case m=2, n=2 which occurs, for example, in 
the study of coupled electric circuits possessing inductance, resistance 


and capacity. In this case each element of 
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cis 
— N(z) 
is a proper rational fraction unless det P(s) is a degree less than 2. Since 
the number of arbitrary constants involved in the solution of the vector 
differential equation is precisely the degree of det P(s), it follows that 
if any element of 
BS ai 

—— N(s) 
is not a proper rational fraction the number of arbitrary constants 
involved in the solution is at least three less than the number (4) that 
would be involved if the vector differential equation were normal. 

Since, then, each element of 


De 

— N(z) 
is, by hypothesis, a proper rational fraction, there exists an mxm matrix 
R= R(t) whose Laplace Transform is 

1 

— N(z) 


- 
~ 


and we have 








_. R(+0)  DR(+0) D*'R( +0) 
NG) -" so = 


or, equivalently, 





DR(+0) D?*R(+0) 
+ 


° oe 
~ - 


N(s)=R(+0)+ 


On multiplying (on the left) this retation by P(s) =Cos"+-+--+ Cy, we 


obtain the series of relations 
( ‘oR ( + Q) = 0; 
C,DR(+0)+ C,R( +0) = 0; 


CoD"R( +0) + «+ +CnR(+0) = En; 
CoD"™*!R( +0) + ---+C,DR(+0) =0 
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In othes words: 
Each column vector of the mxm matrix R= R(t) satisfies at t= +0 
the initial conditions that must be satisfied by a rest solution of the vector 


differential equation. 
It follows that 


; | Ba 
L\ P(D)R} = P(2)L(R) = P(s)—N(z) =— =L(E,,) 





so that P(D)R=E,,. 


In other words 
R=R(t) is the rest solution of the matric differential equation 
P(D)R=En. 
If r=r(t) is a rest solution of the given vector differential equation 
P(\D)y=f we have 
L(r) = N(s) L(f)=L(R) s Lif) =L(R){LD(f)+f(+0) 
and so r=Rf(+0)+(R * Df) 


where we understand by (R * Df) the following 
t 
(R * Df) = R(r)Df(t—r1)dr. 
0 


This is the Principle of Superposition (of rest solutions) for a vector 
differential equation. That r=r(t) defined by the formula 


r=Rf(+0)+(R * Df) 

actually is a rest solution of the vector differential equation 

P(D)y=f 
is an immediate consequence of the facts that P(D)R=E,, and that 
each column vector of R satisties the initial conditions appropriate to a 
rest solution. (We assume that f possesses a piecewise-continuous mth 
order derivative when evaluating D"r. ) 

If we multiply the relation 


DR( +0) D?R( +0) 
N(s) = R(+0)+ : + =. 


“ x 





on the right by P(s)=Coz"+.---+C, we obtain the series of relations 
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R(+0)C,=0; 
DR(+0)Qo+ R(+0)C, = 0: 


D"R(+0)Co+--- +R(+0)C, =F; 
D"*!R( +0) Cot ---+DR(+0)C, =0. 


On multiplying, then, the equation P(D)y=f on the left by Ri +0) it 
follows from the first of the relations above that 
R(+0)}C,D"'+---+C,}y=R(+0)f. 
Hence the initial values of y, Dy,---, D’~!y must satisfy the relation 
R(+0)C,D""!y(+0)+ ---+R(+0)C,y( +0) = R( +0)f(+0). 
If the equation is normal (so that det (y#0) R(+0) is the zero matrix 
and this relation is vacuous. If 7 >1 we obtain an additional relation 
connecting the initial values of y, Dy,---,D""'y by multiplying the 
equation P(D)y=f on the left by DR(+0) and adding to this the 
relation 
R( +0), C,D"+---+C,D}y=R(+0)Df 
(obtained by differentiating with respect to ¢ the relation just given) 
and then using the second of the relations given above. If 7 >2, another 
such relation is obtained by multiplying the equation P(D)y=f on the 
left by D?R(+0) and eliminating D”y by using the first three of the rela- 
tions given above and so on. We suppose in what follows that the 
initial values of y, Dy,---, D"~'y satisfy, when the given vector ditfer- 
ential equation is abnormal, these necessary relations. 

We now introduce the mxm matrix Q=DR each of whose elements 
is the derivative with respect to ¢ of the corresponding clement of &. 
Since 

L} P(\D)Q} =L\D P(D)R} =2 L} P D(R)} —} P D(R)§ 120 
=F,,—E,,=0 
we have P(D)O=0. In other words: 


Each column vector of the mxm matrix Q is a solution of the associated 
homogeneous vector differential equation P(D)u=0. 


It follows that each column vector of the matrices DO, D?Q,. >> is a 
solution of this homogeneous equation Hence each column vector 


of the matrix 
So= D"I100C)+-D" QC, 4+ >>> + QC i-1 


=D"RO+ >: +DRC,..1 
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js a solution of the associated homogeneous vector differential equation 


Furthermore 

So( +0) = D"R(+0) Cot - -- +DR(+0)C,_1=En —R(+0)C, 
and DSo(+0)= —DR(+0)C,; D?S(+0) = —D?R(+0)C, and so on. 
Similarly the mxm matrix 

S,=D"—2OCo+ - >» +OC,2 =D" RCy+ + + + DRC 


is such that each of its column vectors is a solution of the associated 


homogeneous differential equation and 
Si(+0) = —R(+0)C, 1; DS\(+0) =E,, —R(+0)C,—DR(+0)C,-1, 


D2S\(+0) = —DR(+0)C,—D?R(+0)C,_; and so on. Continuing in 
this way we construct the matrices S3,---, ending with S,_;=QCo and 
we set up the vector 


y=r+Soy(+0)+.8, Dy( +0) + ---+5,-.D"-1y(+0) 


where 7 is the rest solution of the equation P\ D)y=f. Since each column 
vector of the various matrices So, . 5,1 is a solution of the homo- 
geneous equation P(D)u=0, y is a solution of the non-homogeneous 
equation P(D)y=f. Since r(+0)=R(+0) f(+0) (by virtue of the 
Principle of Superposition) the value of y at t=+0 is 


2 


R(+0)f(+0)+ }E,,—R(+0)C, }v(+0) 


—R 4-0) C, Dy +-()) —..-—k +0)C)D""'y +0 =Y¥ +). 


(for 


ing this argument we see that the vector 


y=rtSoy( +0) +: --+5,.D"y(+0) 


R(+0)} C,D""!y(+0)+ ---+C,9(+0)} =R(+0)f(+0))  Continu- 


is the solution of the vector differential equation P(D)y=f which as- 
sumes, together with its derivatives up to the (7—1)st inclusive, at 
(= +0 the (properly assigned) values y(+0), Dy(+0),--- D"~'y(+0). 


Note. When the vector differential equation is normal the initial 
values of R and of its derivatives up to the (7—1)st, inclusive, are zero 
while D"R(+0)=Co7. It follows that 

Sol +0) =E£,,, DSo( +0) =0,---, D"~'So(+0) =0; 


similarly S;(+0)=0, DS,(+0) =En,->-, D"'Si(+0) =0; 


ths 1( -+-()) =(),- me tal pets; 1\ +()) = F... 
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Since, when the vector differential equation is normal, R(+0)=0 each 
element of N(z) is a proper rational fraction and 


L(Q)=L(DR) =2L(R) —R( +0) =sL(R) = N(z). 


In this case, then, it is easier to obtain Q directly by analyzing the 
elements of N(s) into simple fractions and then R follows from the 
relation 


R(t)= O(r)dr. 


0 


Since L(r) = N(z) L(f) the Principle of Superposition may be put, fOr a 
normal vector differential equation, in the form 


r=Q*f, 
Example 1. 
(D—1)x—2y=t 


—2x+(D—-l1)y=t 


1 0) 


“Dik Thus the vector differential equation is 


Here n=1 and Cy= 
\ 


normal and x(+0), y(+0) may be assigned arbitrarily. Let us set 
x(+0)=2, y(+0)=4. We have 




















(s-1 —2) 
P(z)=|" 5 | , det P(s) =s?—2s—3 =(z—3)(s4+1) 
\ —2 3-1) 
a _ 
| (s—3)(s+1) (—3)(z+1) 
sothat N(z)=P-(z)=| | 
| 2 s—1 
(g—3)(s+1) (s—3)(s+1) | 
a 1 1 1 | 
ie ‘ > ‘ 9/ ‘ ie | 
| 2(3—3) 2(s+1) 2(s—3) 2(2+1) | 
a | 
gar 1 1 1 | 
| + | 


| (2-3) ~ 2(s+1) 2(s—3) 2(s+1) | 
Since L(Q) = N(z) we have 


1,3t 1 o~t 13 lo~t 
2% + 3e 26° —2 | 


1,3t_1,-t 1,3t41,-¢ | 
3e%'—se 2e*'+5e J 


Q(t) = 
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Since f=v(t,t) (by which we mean that f is the vector whose coordinates 
are t, t) we have Q(t—71)f(r) =v(e3"'’7, e387) so that 


t 


r=Q*f= O(t—7)f(r)dr 


0 
= ((1/9)e*#—(1/3)t—-1/9, (1/9)e**—(1/3)t—1/9). 


Since Co= 2 the desired solution is obtained by adding to 7 the result 
of operating on v(2,4) by Q(t), namely, v(3e?!—e7', 3e8‘+e7-').. Hence 


x =(28/9)e8 —e~'—(1/3)t-1/9; y= (28/9)e8+e-'—(1/3)t—-1/9. 
Example 2. 
(2D —1)x+ (8D —2)y=te' 
(2D+1)x+(8D+2)y=te?! 


This is an abnormal first order vector differential equation since the 
9) / 


oa e | 


5s 7" ‘ 2z—-1 32-2) 
is singular. Since P(s) =}. a 


.s \2s+1 35+ 


det P(s) =2z and so 


we have 


matrix Co= 


/ 




















( \ 
| 3 | 3 1 | 
62. eee ae 
N(s) =P-\(s) = | 
\ ] 1 
| <2 = 
22 2s 
( } 
3 + l 3 Ps 4 
= aa a cae oa 
— Ms) | 2% 2s an 
Hence L(R) =———= | | 
‘ il I l - 
he z < 22? > 23? | 
i 3 3 | 
| —Te et | 
so that R(j= 
— owes | 
= _ = - = | 
} 
ine f 3 3 | 
Hence Q(t) = | 1 1 | Since R(+0) = | 9 ~ 9 | 
aia ane | } 
7 a L =f 1 | 
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a ears —92) me —6| ws 
(= g| wehave R(+0)C,=) 5 4} 80 that the initia} 


values x(+0), y(+0) of x and y must satisfy the relation —3x(+0) 
—6y(+0)=0 or, equivalently, x(+0)+2y(+0)=0. Since Q(t—7)f(7) 
=y(r(e7+e2"), —37 (e7+e?") ) we have 


he | oO cae a ] 
QO *f=v | —t+e'(t—1) +—e*(2t-—1), ——— —e'(t—1) —— e*! (2-1 
| 1 S > 8 


Hence r= R(+0)f+0 * f 
5 ao a 5 x Tee 
= a toare (5t—2) “ae (4/+1), a (3t—1 ren (G+ I 
6 | 
—2 -3)] 
ing v(4x(+0)+6y(+0), —2x(+0)—3y(+0)) to v(t). Taking into ac- 
count the relation connecting x(+0) and y(+0) we find 


Since Q(t)Co= the general solution is obtained by add 


5 l l 
x= y +> e'(St—2) “——" 4{¢-+1)+x(+0); 
5) l ae 
y= —— —— e'(3t-—1) + — e*"(6t +1) +9(4+ 0). 
S Z S 


Note. In this example f(+0)=0 so that, even though R(+0)#0, 
r(+0)=0 since 7(+0)=R(+0) f(+0). Hence y is found by adding 
v(x(+0), 9(+0) tor. 


Example 3. 
(2D?—D+9)x —(D?+D+3)v=0 


(2D?+D+7)x—(D°?—D+5)y=0. 


(9 in 
Here (y=), so that the second order vector differential equa- 
jon is al 1 P(e) = (227 2+9 —(s*-+2+3)) 
tion 1s abnormal. (3) = x es mE os, 
(22°-+-2-+ 4 — (z?—2-+-5) } 


det P(s) =62*—6s°+ 243 —24=6(s—1)(2°+4). Hence 














—(z?—s+5) s°+z2+3 

1 | (s—1)(32-+4) (s—1)(s2+4) | 

N(s)=P"'(2) =—— | 7 , 
| —(2s?+2+7) 23°—s+9 

(g—1)(s?+4) 
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l l i space 
1 | s°+4 s—1 3°+4 s—l | 
ee { 
6 | 2 l 2 | 
). stt+4  3—] gtd peer 


Since each element of N(s) is a proper rational fraction we have 
Ri +0) =0, 





3 ; DR(+0) 
isince N(s) = R(+0)-+—— + :--; , and so L(Q)= Nz). 
\ - / 
/ 
I 
| = sin 2/—e' = sin 2/+e' 
Hence (QU(t)=— 
6 
—— sin 2¢—2e' — — sin 2t+2e' 
| 2 2 
} 
: Ss 1 Il | 
saa haes — — (cos 2t—1) —(e‘'—1) Fos 21) be! 
' j ; 
R(t) =—| 
6 
pe a Dae 1) > (cos 2t—1)+2(e'— 1) 
Since 
, (—1 1) . {-1l 1) . {9 3 
D R(+0)=0(4+0) = ,» C= |, C=), nh 
bie tied [~9 2] oo 1) “(7 —5) 
(2 2) . {-2 —2) 
wehave  D R(+0) C= s)> DR(+0)G=| _| {| > Hence 
- Se ay 


x( +0), v(+0), Dx(+0), Dy(+0) are connected by the relation 
— 2x(+0) —2y(+0)+2Dx(+0)+2Dy(+0) =0 
or, equivalently, 
x(+0)+9(+0) —Dx(+0)—Dy(+0) =0. 


Since f is the zero vector, r(t) is the zero vector. Again 


,( Zun® —sin 22) 

: | 

VOo= = 
6 


[ —2 sin 2 sin 21) 
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1 { cos 2t-e cos 2t+e') 
DQ=— { | 
6 (| —cos 2/—2e' —cos 2/+ 2¢') 
, 1 ( 4 cos 2t —2 cos 2/) ’ l (2e' Je! 
6 \—4 cos 2 2 cos 21) 6 \4e’ de! 
: ’ | ( 2 cos 2t+e' —cos 2t+e! 
so that DOC ot OC, = — | 
3 \—2 cos 2¢+2¢! COs 21-4 2¢') 


l 
Hence x(¢) < 1 (2 cos 2¢-+e')a(+0)+(e'—cos 2¢)y(+0)} 


l , : 
+—}2 sin 2¢ Dx(+0) —sin 2¢ Dy(+0)}; 


6 
‘ t ( 
y(t) orl (2e'—2 cos 2t)x(+0)+(2e'+ cos 2t)¥( +0) ; 
Bg ; By dk i 
+e —2sin 2t D x(+0)+sin 2t D y(+0);. 
) 


5. Ihe determination of functions from their Laplace Transforms. 
It is necessary for the solution of problems by the method of Operationa|! 
Calculus to be able to determine the function f(t) whose Laplace Trans- 
form $(3)=L}f(d)} is in our possession. The following theorem is use- 
ful in this connection: 


Let (zs) be an analytic function and let a be a complex number (whose 
real part is positive) which is such that 2%(s) ts analytic at s=@ so that, 
near = ©, 








l -«< he 
= 
9! a ze on 
~ 0 ~ 
= a : 
Then the power series ¥ ———— converges for every t and the func- 
0 '(n+a) 
ce 6) C = 
: P - n a 
tion f(t)=te' > ———-; ¢>0 
. r(in+a) 


has, at points s whose real parts are sufficiently large, (2) for its Laplace 
Transform (Tp) being the Gamma function which is defined by  th« 
formula 
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@ 


I'(p)= e~1?-'dt, R(p)>0. Thus P'(p+1)=p I\p)). 


it) 
° 


A proof of this theorem may be found in the book referred to in 
footnotel, 
On setting a=} we obtain the following special case of the theorem 


just given: 


\ 











' , l (| Cg 
The function of which - o+-—+ —- +: | is the Laplace Trans- 
form i f(t) sr eee. 
orm is (t) = ——————-¢ ¢ ie 
I= Taye | “T1772 *1/2-3/2 
( cosh # ) r(3) ( T ‘ 
In particular LL) | = eens abet ‘ened ide 
( # 3! ez 
Again, on developing (14 3?)7! near s= @ we obtain 
tiie ——— 
iia (1+3s7)3 


where Jo(t) is the familiar Bessel function of the first kind of zero order. 
Similarly on developing 


] 
3(1+37) 





near 3= © we find that 


Lierf #@) =—_—_———_ 
2( 1+ 37) 


where erf x is the error function defined by 


z 


e~” dt. 


2 
erf x=— 
r 


—\/2 
near s= @ we find that 





On developing 


e 1 





L{ Jo( 2t3) } = 


~ 
r=) 
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and similarly we obtain the more general result 









1/2 
( 
tof 9 ‘ 1] ~ 
Lyte Ja_1(2t)*| =—— ,  R(a)>0. 


oa 
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COLLEGIATE ARTICLES 


Original Papers Whose Reading Does Not Presuppose 
Graduate Training 


A Novel Algorithm at the Freshman Level 


by GERALD B. Hurr 
Southern Methodist University 


Introduction. It is a tedious task to tabulate the values of a poly- 


nomial such as 





(1) 2x4—4x3— 10x7+ 11x+7 


by using just the Remainder Theorem and synthetic division. The 
following algorithm avoids the copying and recopying of the coefficients 
and the determining of unused quotient coefficients. 

Having determned —2 as a lower limit to the zeros of the poly- 
nomial (1), the work indicated below is performed: 
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2 —4 —10 11 7 
D, —4 16 —12 2 
2 —sS 6 —] 9 —2 
—2 10 —16 
2 —10 16 —17 —|] 
0 0 
2 —10 16 0 
3 
2 — I 
The reader will observe that the polynomial was divided by »—(—2 


the quotient by x—(—1), ete., until the operation reached a natural 


conclusion. The work is continued as shown below: 














ye =¢ 
3 ——— —17 
a aa sae 
32 15 
16 15 a Q 
—24 —16 —] 
i =. =f J F 
s 0 —16 —1i 
2 OO ai = 27 -% 2 
8 24 32 15 
2 8 16 a 3 
In this process the rows opposite —2, —1, 0, 1,-.- are computed from 


the top down by multiplying the entry below the coefficient of «* by k 
and adding the result to the entry in the row above and the next column 
to the right. The numbers 9, —8, 7, 6,---, which occupy the positions 
of remainders in synthetic division, are the values of the polynomial at 
—2, —1, 0, 1,---. There is an instantaneous check at x«=0, and a 
check by the usual method at x=3 checks that value and also indicates 
that the intermediate values are correct. 

The writer has found this algorithm useful as bait to be thrown 
to promising undergraduates. It is apt to lead them to ideas not en- 
countered in the usual courses and yet well within their reach. The 
following paragraphs give one explanation of the algorithm and indicate 
its connection with familiar problems. 


1. The construction of the algorithm. In formal work in the ele- 


mentary calculus of finite differences, the polynomial x(« —1)- --(«—n+1) 








140 © 





MATHEMATICS MAGAZINE 













plays an important role and is designated by [x]". The usefulness of 
the symbol [x ]" stems from the relation 







(1.1) A[x ]*=n[x]"", 
where Af(x) =f(x+1)—f(x). To change a given polynomial 
(1.2) f(x) Saux" +ayx"!+---+a,ix+a, 






to the form 


(1.3) 






bo[x—a]"+), [x—a ]"'+ a6 +bD,_1 [x—a ]+2,, 
note that 0, is the remainder when f(x) is divided by x—a@ and the quo- 
tient Gn_1(x) is bo[x—at1 ]"'+),[x-—a4+1]"?+---+,-;. Thus 6,., is 


the remainder when g,_;(x) is divided by x—a+l1 and all the numbers 
bn-2, On_3,- > +, bp may be so determined.* 





















(1.4) The polynomial ayx"+a,x""'+---+a,.x+a, may be put in the 
form bo[x—a }' +b, [x—a J! + --- +0, 1[x]+0, by the algorithm: 




















ay ay ac An-1 ay 
ado ees aqn-2 Ad n-1 
qo Qo Gu. | b, a 
(at1)qo - a (at1)q'n-2 | 





qo. qi’ ee | Dna 


The numerical work (2) in the introduction has the effect of estab- 
lishing the identity 


2x4— 11x?—10x°+11x%+7 
=2[x+2]}'*—8[x+2 P+16[x4+2 ?—17[x+2]+9. 


The completing of the algorithm as carried out in (3) of the introduction 
is now easily explained. 


For f(x) in the form (1.3), Af(x) is given by 
Af(x) =nbo[x—a J" 3+ (n— 1b [x—a J r+ ++ +0,1, 
and f(x+1)=f(x)+Af(x) =bo[x—a ]"+(b:+nbo)[x—a]' 3+ --- 
+ (bn-1 + abn_2) [% —a J+ (bn +bn-1). 
Hence, setting x=x—1, 
(1.5) f(x) =bo[x—a+1 + (b,-+nbo)[x—a+1 P+: -- 
+ (byt 2by-2) [x a+ 1+ (bntbp-). 


*Whittaker and Robinson ‘‘The Calculus of Observations”, Chap. I. 
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This means that if f(x) is in the form (1.3) for a given a, it may be put 
in the same form for a=a+1 by the equations: 


by=bo; b, =b,+nbo; i ae ba-1=dp 1+2bn_2; bn =b,+by-1. 


Inspection of (3) of the introduction now shows that 2x4—4x3—10x? 
| |lx+7 was put successively in the form (1.3) fora=—2, —1,0,1,---. 
The coefficients of these forms must be read diagonally from lower left 
io upper right and the last coefficient, },, must in each case be f(a). 

2. An inverse of the algorithm 1.4. We will say a polynomial is in 
the a-form if it is in the form (1.2) and in the }-form if it is in the 
form (1.3) witha=0. The algorithm (1.4) with a=0 provides a method 
for going from the a-form to the b-form. What about an algorithm 
for going from the 6-form to the a-form? 

To set this up in a simple manner seems to require the following 


lemma: 


2.1) Jf the numbers By, By,- ++, By1, By are given in the identity 


By[x—1 $+ By[x—1 J '+---+B,1[x-1]+B, 
=bo[x }*+),[x I+ - --+),1[4]+5, 


the numbers by,by,: + + ,bx-1,0% may be computed by 





Be yf Bs By 
—kby athe oe —2h,. 9 —b, 1 
bo bh, soe Oy b, 


Setting a=0 and n= in the identity (1.5) yields 
b [x }§¥+), [x i+ ons +b,_i[x ]+, 

=bo[x—1 + (b,4+ Rbo) [x —1 I+ - ss + (bg 1 + 204-2) [x — 1 J+ (b,4+0,-1). 
Thus the identity (2.1) is equivalent to the equations 


Bo=bo; By=bi+kbo;---;  Bes=b,-1+2d,_2; By=by +d, 


bo= Bo; b= Bi-—kbo;-- +5 Oy = Bes —2dy_2; b= By—dy-1. 


The algorithm is merely a schematic solution of this system. 
Now consider the problem of determining the numbers a,, 
@n1,'**, Qy, Ag from the identity 


(2.2) bo[x '+4, [x J+ de ae +b,_1[x ]+0, 
=Aypx"+0,x" 1+ ---+4,_1x+4,. 
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Clearly a,=6,. Cancelling these and dividing both sides by x yields 
bol[x—1 P+), [x —1 JP 2+ - On) Sax" !+ aye" 2+---+a,-. Now 


by the algorithm of (2.1) this may be rewritten 
bo’ [x J’! +,’ [x 2+ + +5", Sage" aye" 2+ + +a. 
Thus a,_1:=0’,_; and there results a new identity 
bo [x—1 J"? +), [x —1 3+ + + +0’, 9 Haye" 2+ ayn" 34+ +. +a, 9, 
which may be transformed by (2.1) into 
bo’ [x PP + by [x [8+ 4b", op Haun" 2+ ax" 34+ +++ +a,-9. 
Thus @,_2=b",_2 and this process may be continued. 


(2.3) The algorithm for going from the b-form to the a-form is : 








| ’ 
bo b oe Dy_2 Dna | b, =a, 
‘ | 
—(n—1)b,y' — 20" 53 —b',_2 
/ / / 
bo by tl b n—2 | ” 1—4yn_-1 
—(n—2)b,”’ —b" »_3 
bo bh,” at, & | bi 2 =n 2] 


or example, [x ]}'+15[x ]'+65[x >+90[x }??+31[x]+1 is coverted into 
x§+ 5x4+ 10x3+ 10x?+ 5x+1 by 














1 15 65 90 nt 4 
1 11 32 2% | 5 

1 8 16 10 

1 6 | 10 

1 5 


Here, and in later numerical examples, the intermediate rows are de- 
leted. 


3. Two applications. The sum of a finite series whose general 
term is a polynomial. 


The formal relation A[x ]'=n[x ]’~! means that a polynomial in the 
b-form may be summed easily. Consider the sum defined for any posi- 
tive integer x: 


z 


(3.1) S(x) =o f(i). 


i=] 


Clearly AS(x)=f(x+1) and f(x+1) may be put in the b-form by 
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From f(x+1) =bo[x "+O, [x J'+ --- +b,-1[% ] +0, 


the algorithm (1.4), 
we conclude that 


S(x) = bofx J08/(n +1) +), [x J! /n+ ++ +b, [x 2/2+0,[x +c. 


Since S(1)=f(1)=),, we conclude that c=0 and 


b ox |" +1 b, [x |’ bo-1 [x P 
2) S(x) = Pe Mea: et 
.2) \ mes" + - y saat 9 





+b, [x ]. 


~~ 


It would be easy in a numerical case to use the algorithm (2.3) to 
put the result in the usual a-form, 

Undergraduates are frequently fascinated by formulas for the sums 
of the power of the integers from 1 to m. A particular case of such a 
problem is worked out below by the method outlined above. 


n 
S;(n) = ry 1° 
i=1 


AS;(1) = n°+5n4+ 1003+ 1002+ 5n+1 
=[n +15[n }4+65[n P+90[n P+31[n” ]+1. 
31[n}? 


» 


_ 





S.(n) = Sad 43[n P+ 430[” P+ +[n] 
) 


65[n ]4 
( 4 


- 
= —— (2[n P+36[n P+195[n }$4+360[n B+ 186[n 2+ 12[7 ]) 



































12 
l . — 
-—< (2n°+6n5+5n4—n?), 
1 5 10 10 5 | 1 
l 6 16 26 | 31 
1 8 32 | 90 
I 11 65 
1 15 
2 36 195 360 186 12 | O 
2 26 87 12 | 0 
2 18 13 | -1 
2 12 | 0 
2 8 | 5 


bo 
for) 
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Interpolation. The algorithm of (2) and (3) in the introduction 
enables one to compute the values of a polynomial f(x) at a set of evenly 
spaced values of x. If on the other hand, the values of a polynomial of 
degree m are given at +1 evenly spaced values, the coefficients of f(x) 
may be computed by reversing the operation. 

The algorithm for changing from the b-form to the a-form fits in 
nicely with the usual theory. If f(x) is represented by the Gregory- 
Newton formula 


A’f(0) A"“1f(0) 
(3.3) f(s) = ert AO) 


[x J’-'+ --- +Af(0) [x ]4+-f(0), 
(n—1)! , : 
where A"f(0) = [A"f(x) ],-0, then the b-form of the polynomial is known, 
once the successive differences at zero are computed from a table of 
values. The a-form then follows easily by the algorithm (2.3). An 
example follows: 

Given the table of values f(x) =3,10,17,33,91 at x=0,1,2,3,4 the 
differences may be computed as below: 


3 











where f(0), Af(O0),--- are 3,7,0,9,24 


24 Ofx 7x 
Thus f(x)=—Z Ee i+ i 4 te) 43 





= 4 (2[x ]*+3[x P+ 14[x]+6) 
= 1(2x4—Qv3+ 13x?+ 8x+6). 


Southern Methodist University. 





CURRENT PAPERS AND BOOKS 


Edited by 
H. V. CRAIG 


This department will present comments on papers previously pub- 
lished in the MATHEMATICS MAGAZINE, lists.of new books, and book re- 
views. 

The purpose and policies of the first division of this department 
Comments on Papers) derive directly from the major objective of the 
MATHEMATICS MAGAZINE which is to encourage research and the produc- 
tion of superior expository articles by providing the means for prompt 
publication. 

In order that errors may be corrected, results extended, and interest- 
ing aspects further illuminated, comments on published papers in all de 
partments are invited. Comments which express conclusions at variance 
with those of the paper under review should be submitted in duplicate. 
One copy will be sent to the author of the original article for rebuttal. 

Communications intended for this department should be addressed to 


H. V. CraiG, Department of Applied Mathematics, 
University of Texas, Austin 12, Texas 
Sequential Analysts of Statistical Data: Applications. Prepared by the Statistica! 
Research Group, Columbia University for the Applied Mathematics Panel, National 
Defense Research Committee, Office of Scientific Research and Development. Colum- 
hia University Press, New York, September, 1945. $6.25. 


This collection of seven stapled pamphlets in a ring binder deals with a new statis- 
tical technique, that of sequential analysis, introduced in 1943 by A. Wald for use in 
war research. An earlier ‘‘Restricted’’ edition for limited circulation preceded this 
printing and has, no doubt, served to reduce errors and misprints common to photo- 
vraphed typewritten copy. 

No attempt is made to present details of the mathematical theory. Indeed, in 
the appendix which is mathematical in nature the following statement occurs: ‘‘This 
ippendix covers only a small part of the theory covered in SRG.75, and that part is 
covered briefly andintuitively. Any mathematician who may stray into this Appendix 
should be assured that the validity of the conclusions in no case depends upon the 
type of reasoning presented here.’’ By way of explanation, SRG 75 is Statistical 
Research Group report number 75, Sequential Analysts of Statistical Data: Theory by 
\. Wald, September, 1943. A more accessible paper for the general reader is Wald’s 
paper, ‘‘Sequestial Tests of Statistical Hypotheses,’’ Annals of Mathematical Statistics 
Vol. 16 (1945) part 2. The mathematician is therefore urged to read tis paper and 
similar research articles for proofs. 

The style is such that users who are neither experts in mathematics nor statistics 
can apply the methods developed to any one of the specific applications described. 
Thus, the second pamphlet deals with a binomial distribution. Suppose, for example, 
that a group or lot of a particular type product is being examined to see if it is a good 
or a bad lot. In many cases, the test might be destructive to an individual item, or 
even if non-destructive too costly to apply to each individual item. Prior to the in- 


vention of sequential analysis it was customary to select a sample of size N, test these 
V items, and then either accept or reject the lot from which the sample was drawn. 
In sequential analysis, the following procedure is used: After each observation, 


one of the following three steps is taken: 
(1) the lot is accepted, or 
(2) the lot is rejected, or 
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(3) judgment is deferred and another observation is made. The nature of the 
theory is such that good lots are quickly accepted, i. e. alternative (1) can be accepted 
before N trials which would have been necessary under the old procedure, bad lots are 
quickly rejected, but border-line lots still require extensive testing, perhaps more 
than N items. 

If the probability of rejecting a lot with p; proportion of bad items is a and if the 
probability of accepting a lot with p, proportion of bad items is 8, then Wald found 
that the criteria for choices (1) and (2) above are 


pi! 1—p, n—d Ea 





(1) for acceptance 
d o-4 = ; ia P 
"(1 — ps) p 
a ; y 
(2) < —— for rejectance, 
l —p 


where m is the number of items tested to date and d is the number of bad items. Test- 
ing continues until one of these relations is satisfied. 

Both numerical and graphical procedures useful in keeping track of the progress 
of the test are explained, and tabular data and nomographic charts with full explana 
tions for their use are given. 

The removal of the ‘‘Restricted”’ tag from this statistical problem is a welcome 
move and although the volume under review is essentially non-mathematical, the 
related mathematical portions will enrich the literature of mathematical statistics 


The University of Texas. ROBERT E. GREENWOOD 


The Magic of Numbers. E. T. Bell. New York: Whittlesey House, McGraw- 
Hill, 1946. X+418. (Frontispiece). $3.50. 


While considering Pythagoras’ numerology Bell repeatedly asks the question 
were numbers discovered or invented? The book doesn’t reach a definitive answer t 
the question but the case for discovery seems more plausible after the reader has plowed 
through the story and met on the way such variant personalities asThales, the fabulous 
Empedocles, Pythagoras and his Brotherhood, Plato, Dante, Bacon, Galileo, Newton, 
Berkeley, Saccheri, Lobachewsky and Kant. Integrated all this forms a background 
for the controversy generated by Sir Arthur Eddington from his claim that the funda 
mental physical laws can be deduced in an epistemological manner without recourse to 
observation or experiment. This controversy receives only limited space in the last 
chapter but deserves much more. Being in the foundations of physical science it is 
of transcendent importance in this atomic age. 

A highlight of the book is the passage where Bell considers an infant subjectivels 
learning Boolean algebra in the relation J and not —J. The book is written so only 
ordinary arithmetic is prerequisite to its reading, the mathematics of the book being 
limited to the beginnings of the theory of numbers while touching upon formal logic 
and non-euclidean geometry. 

For reference while reading the book it is well to have Runes’ Dictionary of Philoso- 
phy on hand. Helpful background information on almost everything from Absolute 
to Zeno is thus available and Alonzo Church's formal logic definitions are outstanding 

For more information on Eddington’s ideas on epistemology one should consult 
The Philosophy of Physical Science, 1939. 





Colorado Springs, Colorado. DUANE STUDLEY. 





HISTORY AND HUMANISM 


Edited by 
(>. WALDO DUNNINGTON and A. W. RICHESON 


Papers on the history of Mathematics per se, the part it has plaved 
in the development of our present civilization and its relation to other 
sciences and professions are desirous for this department. 


A French Mathematician of the Sixteenth 


Century 
Jacques Peletier (1517-1582). 


bv V. THEBAULT 


Tt nnte Sartne France 


Biography. Jacques Peletier, French writer, poet and mathemati- 
cian, was born at Le Mans, July 25, 1517. He was sent at an early age 
to Paris and there he attended Navarre college, where bis brother Jean 
taught philosophy. Soon after leaving college he prepared to become a 
teacher. In 1544 he obtained through the influence of the bishop of 
Le Mans,René du Bellay, the position of tutor of mathematics at the 
college of Bayeux. Three vears later, in 1547, he was at the head of 
this school, when he delivered from the puplit of Notre Dame the funeral 

ration of Henry VIET of England. This event seems to have been a 
turning point in his life. 

The restlessness of Peletier prompted him. to resign, in order to go 
to live in the house of the printer, Vascosan, where he conceived the 
curious project of reforming spelling according to pronunciation.* This 
forms the subject matter of his Dialogue on French Spelling and Pro- 
nunciation, in two volumes. 

Peletier had scarcely completed this work when he became deeply 
interested in medicine, which he went to study in Poitiers in 1550. He 
staved successively in Bordeaux, Lyon and Rome. He then seemed to 
be tired of wandering, and settled down a while in Paris, where he re- 
ceived his Licentiate in medicine. The turmoil of the Wars of Religion 
again took him away from Paris. After a long stay at Annecy, Savoy, 
he returned to Paris in 1578 to become the head of the college of Le Mans. 

Since he had always remained faithful to the Church of Rome, in 
an era when religious differences menanced French national unity, the 





Faith that he had inherited from family tradition caused him, toward 
the end of his life, to humbly implore divine aid in contemplating the 


*Dialogue de l Orthografe e Prononctation francoese. \tarnef, Poitiers, 1550. 
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eternal laws which govern numbers and figures. It is in his beloved 
geometric language that he phrased his last prayer: 

Thou, One, All Infinite, Thou the Circle and the Center, 

Whence comes every line, where every line returns, 

Make me end in Thee, I, who in Thee began. 


Peletier died in Paris in July, 1582. 


The Mathematical Work of Peletier. To understand the meaning of 
the first mathematical work published by Peletier we must take into 
account the extreme poverty of Western thought in the domain of the 
exact sciences up to the sixteenth century. After the fall of the Western 
Empire all that remained of the science of Euclid, Archimedes and ap- 
pollonius consisted of meager summaries, simple enumerations of the 
most elementary rules transmitted by Boéce and Cassiodore. Arithmetic 
and geometry were reduced to statements covering a few applications. 
Writers on the subject of algebra, which had come from India by way 
of the Arabs, were reluctant to solve problems in a general manner; 
they repeated the reasoning for each particular problem. Furthermore, 
no symbolism had been adopted which permitted a common treatment 
of all numbers, positive or negative, integers or fractions, real or complex. 

From 1515 on, the wars in Italy had brought French scientists into 
contact with the students of the learned Greeks who, chased from 
Byzantium sixty years before by the Turks, had fled to Sicily and had 
gone little by little toward the North, bringing with them the texts of 
the great mathematicians of the Alexandrian period. It is difficult to 
depict the astonishment of Europeans upon reading such majestic and 
clear monuments of pure logic as are, for instance, the Elements of 
Euclid, the Conics of Appollonius and the 7reatise on Spirals of Archi- 
medes. Thus Peletier, who doubtless had never known of the existence 
of such treatises before his arrival in Paris, gave himself no respite until 
he had learned Greek by himself (Greek was not taught in the schools of 
France until 1532) and had read the original texts. 

For the study and teaching of arithmetic he had at his disposal 
only the manual published several years before by Gemme of Friesland, 
professor at the University of Louvain. He corrected many of the errors 
which abounded in the original edition, and, encouraged by his students, 
published in 1545 a first Latin edition, Its success was considerable, for, 
a rare event for this period, Peletier’s Arithmetic was reprinted, withou! 
change, but in French, at Poitiers in 1548 and 1552, at Lyon in 1554 and 
1570, in Latin in 1568, and finally in Italian in 1567. 

During this time Peletier undertook the study of algebra. Proba- 
bly he was the first in France to publish a treatise in algebra, in 1554 
inaugurating the use of literal symbols, which were later made popularly 
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Viete. He showed that quadratic equations admit not one, but two 
roots, Whose product equals the constant term. This was extended later 
by Viete to cubic equations and by Harriot to equations of arbitrary 
degree. Also he carried over to biquadratic equations what he had 
proved for quadratic equations. 

By publishing his Algebra Peletier became a leader and a popu- 
larizer, but in spite of the clarity of his exposition he was unable to over- 
come the hostility brought upon him by his use of the French language. 
But he was of peaceful temperament, and the only revenge he took was 
to translate the second edition of his Algebra (1560) into Latin, with the 
ironical title, De occult a parte numerorum quam Algebram vocant. 

After 1554 Peletier, without neglecting the new editions of his 
works, devoted himself more and more to consolidating the theory. At 
the same time his enthusiasm for Greek mathematics began to wane 
\fter Cardan had called his attention to flaws in the majestic edifice 
of Euclid’s Elements, he examined closely the foundations of geometry, 
and believed, quite some time before Lobachewsky and Riemann, that 
he had discovered defects. 

Peletier reduced to two the paralogisms and contradictions which 
he believed he had found in Euclid’s Elements, at that time the holy 
hook of mathematicians. The first referred to the proof of the congru- 
ence of figures by superposition and the second to the magnitude of the 
angle formed by a curved line and its tangent. For the first question he 
let himself be convinced by Jean Buthéon or Borel. But for the second 


question he maintained his position. 


In the thirteenth proposition of the third book of his Elements: 
luchid considers the figure formed by a circle and one of its tangent 
lines. By analogy with the angle formed by two intersecting straight 
lines, he considers the angle of contingence and shows that it is smaller 
than any assignable angle. 

Peletier states: ‘‘When I examined this proposition more care- 
fully, it suddenly appeared to me that geometry could not be satisfied 
with it without admitting inconsistency. First, we can not suppose that 
there can exist a quantity smaller than all others, as this angle of con- 
lingence or contact would be, for quantity is what is formed of parts, 
and it is only with respect to quantity that things are said to be equal or 
unequal. On the other hand, continuous quantities are indefinitely 
divisible. This is why it seems to me impossible to reconcile what is 
said here about the angle of contingence with the first proposition of the 
tenth book.”’ 

However, the disturbance stirred up by these discussions as to the 
angle of contingence kept increasing. Clavius of Bamberg, a noted 
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astronomer, had just published in Rome a treatise on geometry in which 
he described the ideas of Peletier as sophisms and hallucinations. Pele- 
tier’s answer only served to spread controversy, as the whole Jesuit 
order backed Clavius, himself a Jesuit. It was to last long after Peletier’s 
death and maybe it would still be going on if Newton himself had not 
definitely settled the argument in favor of Peletier by calling attention 
to the possibility of variation of curvature near the point of contact. 

This controversy caused Peletiec, just a few months before his 
death, to write a short 14-page memoir entitled De contactu. Linearum, 
of which only one copy remains. But this work, by the precision of its 
definitions, the clarity of its reasoning, its careful and well ordered ex- 
position, belonged already to the Cartesian era, fifty years ahead of 
contemporaneous mathematicians. 

Peletier, who had understood perfectly the core of the difficulty, 
analyzed at length the notion of angle. He called attention first to the 
fact that since geometry is the science of magnitudes it can have no 
other object than the study of quantity. Since angles properly belong 
to it, angles must be measurable quantities. Their measure can be 
sought neither in the length of the sides, which can be extended without 
changing the angle nor in the area included between the sides, since 
this area is either unbounded or is that of a triangle if a suitable third 
line is drawn. The magnitude of an angle must exist at the vertex 
itself and it will be smaller the more nearly the sides coincide. 

Going from this to the notion of contact Peletier notes that the 
angles formed by homologous sides of two similar figures are equal. It 
follows that angles of contingence, at least those that circles of different 
radii form with a common tangent line, can not have different magni- 
tudes. Therefore angles of contact are not quantities and cannot, 
whatever Euclid may have said, be compared to angles formed by 
inte.’secting lines. 

In conclusion let us cite the following lines written by Abbé M. 
Thureau in his book entitled Jacques Peletier, Sixteenth Century Mathe- 
matician of Le Mans, 1934: ‘In all the abundant scientific literature of 
the Rennaissance I do not know of a single work which even nearly 
approximates Peletier’s De contactu. linearum, a masterpiece incon- 
testably worthy of appearing in an anthology of mathematics if, as | 
hope, one is ever published.” 

Note: In the time of provincial France the city of Le Mans, the capital of the 
province of Maine situated in the west of France, was the birthplace of another illus 
trious mathematician, the Reverend Marin Mersenne, born September 8, 1588 in a 
town only a little distance from Le Mans. He died in Paris, in the convent of the 
Minimes, on rue Royale, September 1, 1648, in the presence of his friend Gassendi. He 


was permanent secretary of the Academy of Sciences at its beginning. His life and 
works are well known. 


Translated by Miss Marie Edith Byrne. 








TEACHING OF MATHEMATICS 


Edited by 
Josep Srrpiin, L. J. ADAMS and C. N. SHUSTER 


This department is devoted to the teaching of mathematics. Thus 
articles on methodology, exposition, curriculum, tests and measurements, 
and any other topic related to teaching, are invited. Papers on any sub- 
ject in which you, as a feacher, are interested, or questions which you 
would like others to discuss, should be sent to JosEpH SEIDLIN, Alfred Uni- 
versity, Alfred, New York. 


Proofs of the Fundamental Theorems of 
spherical Trigonometry 


by C. E. CLARK 


This note presents derivations of some fundamental formulas for 
the spherical triangle. The derivations are shorter and simpler than 
those found in the textbooks. Furthermore, the first proof enables a 
teacher to present the essential applications of spherical trigonometry 
after less than one lecture on the spherical triangle. The formulas 
derived are the law of cosines for sides, the law of sines, the law. of 
cosines for angles, and Napier’s rules. The derivations are shorter 
and simpler than those given in the textbooks for the following reasons, 
The use of solid geometry including the theory of the polar triangle is 
avoided. The only formulas from plane trigonometry used are the law 
of cosines, the reciprocal relations, and the Pythagorean relations. In 
no case is a proof divided into two or more cases depending upon the 
quadrants of the angles involved. Another important feature of the 
discussion is that it presents a brief course in spherical trigonometry. 
This brief course arises from the fact that the first derivation is for the 
law of cosines for sides. This law alone will handle applications of 
spherical trigonometry to the terrestrial and celestrial spheres. 

To derive the law of cosines for sides we denote the sides of a spheri- 
cal triangle as well as their measures by a, 6, and c. The 
opposite vertices and their measures are A, B, and C, respectively. 
The center of the sphere is O. A plane is constructed perpendicular 
to OA. This plane meets OA, OB, and OC in the distinct points, P, Q, 
and R, respectively. Application of the law of cosines for plane tri- 
angles to the triangles POR and OQR gives the equation 


PQ)?+-(PR)2—2(PQ)(PR)cos A =(OQ)?2+(OR)?—2(0Q)(OR)cos a. 
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If no side or angle of the spherical triangie has the measure 90°, the 
lengths PQ, PR, OQ, and OR are replaced by OP tan c, OP tan b, OP sec c, 
and OP sec b, respectively. A simple reduction gives the law of cosines 
in the form 


(1) cos a=cos 6 cos c+sin b sinc cos A, 


That this relation still holds when a measurement is 90° is clear on the 
grounds of continuity. This law of cosines gives three independent 
equations which will solve any spherical triangie. In some cases the 
solution would be long, and the elimination of extraneous solutions 
would be tedious. However, when three sides or two sides and the 
included angie are given, another measurement can be found readily 
by the law of cosines. These cases cover several of the significant appli- 
cations. Examples of such applications are the problem of finding 
distance and bearing when the latitude and longitude of two points on 
the earth’s surface are given, and the problem of finding the time when 
the latitude of an observer, altitude of the sun, and declination of the 
sun are given. 

We turn next to the law of sines. A geometric derivation resem- 
bling the above derivation of the law of cosines is given in the Encyclo- 
paedia Britannica, 14th edition, volume 22, p. 474. This encyclopaedia 
article uses a few results from solid geometry and is not quite as simple 
as the above derivation of the law of cosines. However, it may be 
preferred to the analytic derivation which we now present. For refer- 
ence, we write the following two statements of the law of cosines: 


(2) cos b=cosc cosa+sinc sina cos B, 
(3) cos c=cos a cos b+sina sin b cos C. 
Multiplication of both sides of (1) and (2) by cos a and cos b, respective- 
ly, and subtraction gives the equation 
(4)  cos’a—cos’)=(sin b cos a cos A—sin a cos b cos B)sin c. 
Similarly, multiplication of (1) and (2) by sina cos B and sin } cos A, 
respectively, leads to the equation 
(5) sina cosacos B—sinbcosbcos A 

=(sina cos b) cos B—cos a sin b cos A )cos ¢. 


Squaring both sides of (4) and (5), adding, and cancelling one term 
from each side of the resulting equation, we get the equation 

b 9 6 bd 9 D 
(cos*a —cos’b)?+sin2a cos’a cos?B+sin’b cos’) cos?A 


= sinh cos2a cos?A +sin2a cos’) cos?B, 
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This equation can be written 
(cos2a —cos?h)?+ sin’a cos?B(cos*a — cos?) ) 
— sin?) cos?A (cos*a—cos*h) = 0. 
Removing the common factor and changing cosines into sines We obtain 
ihe equation 
sin’a sin? B=sin2) sin?A, 

Since ali angles involved are between 0° and 180°, we obtain the law 
of sines. 

The law of cosines for angles can be obtained from the theory of 


the polar triangle. However, the geometry involved can be avoided 
by the following analytic proof. The law in question is that 


6) cos C=cos A cos B+sin A sin B cos c. 


To show that this equation 1s an identity, the law of sines is used to 
show that 


: snasin b sin?C 
sin A sin B= 





sin’¢ 


sina sin } 
= — (1—cos?C). 
sin’¢ 





This last expression is substituted in (6). In the resulting equation 
cos C, cos A, and cos B are replaced by the expressions obtained from 
3), (1), and (2), respectively. Then obvious simplifications complete 
the proof that (6) is an identity. 

We can now obtain Napier’s rules from the above results. From 
6) we see that 


7) cos A = —cos B cos C+sin B sin C cos a, 
&) cos B= —cos C cos A+sin C sin A cos 0, 


Henceforth we assume that C=90°. From (3) we see that cos c=cos a 
cos b. From (7) and (8) we get the relations cos A =sin B cos a and 
cos B=sin A cos 6. From (6) we obtain the equation cos c=cot A 
cot B. The law of sines gives the relations sin a=sinc sin A and sin } 
=sinc sin B. Using the relations already derived we see that 

sin a 


tana= 
cos a@ 
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sinc sin A sin B 





cos A 
=sinbtan A. 
From symmetry we have that tan b=sina tan B, Furthermore, 


sin a 
tfana= 
cos a 


tan bcos b 





tan Bcosc 


sin } 





tan Bcosc 


sinc sin B 





cosctan B 
=tanccos B, 


Symmetry gives the equation tanb=tanccos A. Thus Napier's 
ten rules are derived. 

Napier’s rules can be derived from the law of cosines for sides and 
the law of cosines for angles without use of the law of sines. Hence 
Napter’s rules can be obtained by first proving the law of cosines for 
sides, next obtaining the law of cosines for angles by the theory of the 
polar triangle, and finally deriving the rules in the following way. Four 
of the ruses ate obtained as the first four above. Using these 
four results we have 


cos a—cos § cos c¢ 





cos = : : 
sin b sinc 


cOs @—CcOs @ cos2h 





sin b sin ¢ 


cos a sin’) 





sin 6 sin ¢ 


cos @ sin b 





sin ¢ 


cos c sin b 





cos b sinc 


=cot c tan 0. 




















PROOFS OF THE FUNDAMENTAL THEOREMS 155 
This gives the rule tan b=tan ¢ cos A, and using symmetry tan a=tan ¢ 
cos B. Next we have 


sind@=tan da cos a 


II 


tan c cos B cos A 
sin B 





=cot Btan bd. 


This is the rule tan b=sinatan B. From symmetry we have that 
iana=sinbtan A. The remaining two rules are obtained from the 
equalities 


sina=tan hb cot B 


tan bcosc¢ 





cot A 
=tanccosA cosctan dA 


sin A sin ¢. 


Symmetry gives sin b=sin B sin ¢ 


Emory University. 











PROBLEMS AND QUESTIONS 


Edited by 
C. G. JAEGER and H. J. HAMILTON 


This department will submit to its readers, for solution, problems 
which seem to be new, and subject-matter questions of all sorts for readers 
to answer or discuss, questions that may arise in study, research or in 
extra-academic applications. 

Contributions will be published with or without the proposer’s signa- 
ture, according to the author's instructions. 

Although no solutions or answers will normally be published with 
the offerings, they should be sent to the editors when known. 

Send all proposals for this department to the Department of Mathe- 
matics, Pomona, College Claremont, California. 


SOLUTIONS 


5S ) 


Following is a solution to Problem No. 3 (Vol. XXI, No. 1, p. 
proposed by Nev. R. Mind. 


Solution by Laurence A. Ringenberg, Charleston, Il. 


To prove: If the medians of a triangle are proportional to the 
corresponding sides, the triangle is equilateral. 











C ’ &. 
_ ™ z 


£ 
Let AA’=a’, BB’=b',CC’=c’: 


Givena:b:c=a’':b’:c'’. It follows that 





a b oC 2a’ 2h’ 2c’ _a’ Ob c’ 
“Se Bak Gat See” SS i ee ee 


Hence AAOB’'~ABOA' and a=a’. Similarly B=6’ and y=y’. It fol- 
lows that a+f6+y=7/2; a+-B+7'+ABO=r: AB'O=7n/2. Similarly, 
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BA'O and BC’O are right angles. It follows that AA B’B~ACB’B and 
that <A=<C. Similarly <C=<B. Thus the triangle is equi- 


aneular and hence also equilateral. 


Solution for Problem 4, proposed by Pedro A. Piza, San Juan, 
Puerto Rico: 

Solution by Paul A. Clement, University of California at Los 
Anveles. 

Multiplying a'+b2=c? by a?—b’, we get at—bt=c*?a?—c*b?, which 
cives the desired relation a!+c*b?=6'+c?a*. Furthermore, c?b?+a‘ 
c'—c’a®+a'=c1a2%b?, and it can be noted that these relations are 


valid for the sides of any right triangle. 
Taking c=a+l, we find that each function considered above, when 
augirented by 3 is factorable to an algebraic square: 
a*c?+-614+3 =a°c?+(c?—a’)?+3 

=ar°c?+(c+a)?+3 

=a*c?+ 2ac+1+c?+a’+2 

= (ac+1)?+2a?+2a+2+1 

= (a’+al1)?+2(a*+a+1)+1 

—(a’?+a+2)? 
Since 6 and ¢ are integers when a is, we have that a°c?+6!+3 is a perfect 
square when c=a+l. 


Also solved by H. E. Bowie, Indian Orchard, Mass.; Laurence A. 
Ringenberg, Charleston, Ill.; Thomas Griselle, Hollywood, Cal. 


No. 600 (N. M. M.). Proposed by N. A. Court, University of Okla- 
homa, 

Three concurrent non-coplanar lines, bearing the equal segments 
SA, SB, SC, are cut in the points A’, B’, C’ by a variable plane PA’ B'C’ 
passing through a fixed point P. Find the locus of the points common 
to the three variable spheres having A’, B’, C’ for centers and AA’, 
BB’, CC’, for radii 

I Solution by Howard Eves, College of Puget Sound, Tacoma, 
Washington 

Let u, v, w, r, s, t, p be position vectors of A, B, C, A’, B’, C’, P 
respectively with respect to S as origin. Then we have 


U-u=V-V=W'W=R, Say, 
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where k is a scalar. Let x be the position vector of a point common to 
to the three spheres A‘(A’A), BY B’B), C’(.C’C). Then, by the condi- 
tions of the problem, we must have 


(1) lp—r, p—s, p—t| =0, 
(2) x°x—2x-r=k—2u-r, 
(3) x-x—2x-s=k—2z-s, 
(4) x-x—2x-t=k—22w-t, 
(5) r=au, s=bv, t=cw, 


where a, 0, c are scalars Substituting (5) in (2), (3), (4) and solving 
for a, b, c we find 
x-x—k x-x—k x-x—k 


(6) i —— , b= ; 
2(x-u—k) 2(x-v—k) 











Substituting (5) and (6) in (1) we find, after expanding and then remov- 
ing the non-vanishing factor (x-x—k)?, 


(7) (x-x}u,v,w| —2xv-[u| p,o,w| +0] u,p,w| +w|u,2,p] | 
+k[ 2| p,2,2| +2|u,p,w|+2|u,v,p| —| u,v, | =() 


Thus the required locus is seen to be a sphere whose center and radius 
may readily be found from its equation (7). 

The analagous problem in euclidean #-space may be similarly 
solved, yielding an v-sphere for the required locus 


II Solution by Nev. R. Mind. 


The tangent planes a,6,y to the three variable spheres (A’), (B’), 
(C’) at the fixed points A, B, C are also tangent, at those points, to the 
fixed spheres (.S) having S for center and SA =SB=SC for radius, i. e., 
the three planes are the radical planes of the pairs of spheres (.S), (A)’; 
(.S), (B’); GS), (C’), respectively. Thus the fixed point L=a6y is the 
radical center of these four spheres, and their orthogonal sphere (L) is 
the fixed sphere having L for center and orthogonal to the fixed sphere (.S). 

The radical axis of the three spheres (4’), (B’), (C’) passes through 
the point Z and meets the three spheres in the same two points M, M’, 
inverse with respect to the sphere (ZL). On the other hand, the points 
M, M’' ae symmetrical with respect to the plane of centers A’B’C’ of 
the three spheres, these points are therefore equidistant from the fixed 
point P through which the plane A’B’C’ passes, by assumption. Thus 
M, M' lie on the sphere (P) having P for center and PM=PM' for 
radius, and the sphere (P) is orthogonal to (LL). Now the fixed point ? 

















PROBLEMS AND QUESTIONS 159 


is the center of only one sphere orthogonal to the fixed sphere (L), 
hence (P) is the locus of the points M, M’. 

BIBLIOGRAPHICAL Nore. The corresponding question in the plane 
was considered in the Nouvelles Annales de Mathématiques, series 2, 
Vol. 17, 1878, p. 523, Q. 1286.—N. A. C. 


PROPOSALS 


13. Proposed by Julius Sumner Miller, Dillard University, New Or- 


leans, La. 


A uniform chain of length Z and total mass M is held vertically 
with its lower end just touching a platform balance. The fixed upper 
end is released and the chain “accumulates” on the scale pan. What is 
the maximum reading of the balance? 


14. Proposed by A. K. Waltz, Clarkson College, Potsdam, N. Y. 


The tangent of an angle, whose sides slide on two given circles, 
has constant magnitude. Determine and discuss the locus of its vertex. 


15. Proposed by Victor Thébault, Tennie, Sarthe, France. 


A rectangular segement EF of constant length slides on the 
diameter AB=2R of a semicircle (O). The perpendiculars to AB at E 
and F meet (O) at Mand N. Show that the points of intersection P and 
P’of the circles passing respectively through M and N and having their 
centers at fixed points C and D on AB are on a circle concentric to (QO), 
and that, if CD=R, we have the relation 


cos’a+cos*8+cos’y = 1 
between the angles a,8,y made by AM, BN, AQ with AB, Q being the 
point where the circle of center C and radius CA intersects the locus of P. 
16. Proposed by H. E. Bowie, American International College. 


A circle of radius 3-in. is tangent externally to a rectangle at the 
mid-point of one end. Two other circles, both of radius 2-in., are tan- 
gent externally to the two sides of the rectangle and to the first circle. 
The rectangle is 2-in. wide. Find the radius of a circle to which the 
three given circles are tangent internally 


17. Proposed by Pedro A. Piza, San Juan, P. R. 


Let x, y and g be any three arbitrary non-zero integers, and let p be 





160 MATHEMATICS MAGAZINE 


a prime number greater than 2. Prove that for any value of p what- 
soever 
[ (x+y) (x+s) ]?—[x(x+y+z) ]?-[ys]?= peys(xty) (ats) (xt y+s)K 


where K is always an integer. 


18. Proposed by Anonymous. 


Shortage in button-board has led many to substitute sheet rock, 
which comes in sheets 4 feet by 8 feet, in which they have bored holes 
one-half inch in diameter in such a way that no pont on the surface of 
the plaster-board is more than three inches from the perimeter of the 
nearest hole. What is the minimum number of holes that can be used? 


Statin, 
MATHEMATICAL MISCELLANY 


Edited by 
MARIAN E. STARK 


Let us know (briefly) of unusual and successful programs put on by 
your Mathematics Club, of new uses of mathematics, of famous problems 
solved, and so on. Brief letters concerning the MATHEMATICS MAGAZINE 
or concerning other ‘‘matters mathematically” will be welcome. Address: 
MarRIAN E. Stark, Wellesley College, Wellesley 81, Mass. 


Any teacher of mathematics in a secondary school or college who 
has taught five years or less is urged to send this department a state- 
ment of not more than one short paragraph, containing what he con- 
siders the most important advice to give a teacher of our subject as he 
starts his career. The writer’s name should be given, together with 
the name of the institution where he teaches and the number of years 
he has taught. The most interesting statements received will appear 
in this department from time to time. 


One report on a hobby of a mathematics teacher after ‘‘age 65” is 
at hand. A Harry Wheeler of Worchester, Mass., formerly of Clark 


University, has taken up the invention of games, more or less mathe- 


matical. 


Advice needed! What can be done with a college student who 
can do A work but is perfectly satisfied with C ? Such students are 
very infrequent, but exceedingly irritating. Of course, they have a 
right to choose how they will spend their time and energy. But how 
may One inspire them to a better choice than they are making? 
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“The Duodecimal Society of America is engaged in educating the 
public in counting by twelves, instead of tens, and in the application of 
the 12-Base in mathematics, weights and measures, and other branches 
of pure and applied science. Its Annual Awards are conferred on those 
who have made outstanding contributions to the development of duo- 
decimals and to the advancement of the purposes of the Society.” 


New York. THE DUODECIMAL SOCIETY OF AMERICA. 


“The calculation of to many places of decimals has had a fas- 
cination for many computers and various formulae have been used for 
the purpose. In 1841 William Rutherford, using a formula of Euler 
(1764). 

r/4=4 tan7! 1/5—tan7! 1/70+tan“! 1/99, 


calculate the value of w correct to 152 places of decimals. Using the 
formula 
w/4=tan-! 1/2+tan— 1/5+tan7! 1/8 


Zaxharias Dase published in 1844 the correct value to 200 places of 
decimals. In 1853, with the formula of Machin (1706). 


n/4=4 tan“! 1/5—tan™ 1/239, 


Rutherford extended his computation of the value of z correct to 440 
places of decimals. In the same year William Shanks, using Machin’s 
formula, published a value to 607 places of decimals which was extended 
to 707 in a value published in 1873. In 1854 Richter’s independent 
calculation of the value of m showed agreement with the result of Shanks, 
through the five hundredth decimal 

“Using the formula 


w/4=3 tan7! 1/4+tan7! 1/20+tan7! 1/1985. 


Mr. D. F. Ferguson, of the University of Manchester, calculated 2 to 
710 places of decimals and showed (1946) that Shanks’ value for 7 was 
quite incorrect from the 529th decimal place on. Using the Machin 
formula Dr. John W. Wrench, Jr., of Washington, D. C., and Mrs. 
Levi B. Smith of Talbolton, Georgia, recalculated the value of x from the 
500th to the 808th place of decimals. This result which is published in 
Mathematical Tables and Other Aids to Computation, for April, 1947, 
may be regarded as a companion to the value of e, computed to 808 
places of decimals by Peder Pedersen, also published in the above-men- 
tioned periodical, April, 1946. 

“To Euler (1707-1783) is due an extraordinary formula connecting 
the ‘five most important numbers in mathematics’, namely: é74+1=0,. 
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Our notations e, i= ¥—1, are also due to Euler. The use of x as the 
ratio of the circumference of a circle to its diameter was first given by 
William Jones, in his Synopsis Palmariorum Matheseos, London, 1706, 
p. 243.” 


Correction. In the September-October number of the MAGazinr 
on page 59 read Peder Pedersen for Peder Pederson, 


At Duke University an Institute for Teachers of Mathematics was 
held, August 5-15, 1947, with the general theme ‘‘Mathematics at Work”, 
Any reader who attended the Institute is invited to report to this de- 
partment briefly on interesting results there. One or more reports 
could be published. 


Princeton University has published a pamphlet concerning the 
Conference on the Problems of Mathematics held in connection with its 
Bicentennial Celebration in 1946. Very interesting reading! Professor 
Solomon Lefschetz was the director of the Conference. 


A student, desiring to perform the integration 


d® 
tan 9 sec’?@—tan 0 





(1) 


notes that the cotangent is the reciprocal of the tangent, and hence any 
cofunction is the reciprocal of the function. He rewrites the integrand 
as follows 


1 
(2) =cot 0 csc*A—cot 0. 


tan 0 sec’°8@—tan O 





He then performs correctly the revised integration; and produces the 
correct answer. Irritating fact: (2) is an identity. 
Dartmouth College. B. H. Brown. 


Colonel W. E. Byrne of the MAGAZINE'’s editorial staff is going to 
report to this department from time to time “‘news of the present status 
of mathematicians in France, their needs, etc.’’ We learn from him 
that M. Gaston Julia has just returned from lecturing in Sweden and 
Norway. Can any reader supply news of our colleagues in other foreign 
countries? 


Correct Results by an Unorthodox Method. While resolving ‘m- 
proper fractions into an integral function and two or more partial 
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fractions in order to integrate, a student found that after obtaining the 
integral rational part of the answer he could sometimes find the partial 
fractions by an unorthodox procedure. This was to set the original 
fraction 

f(A) A, A, 


oe ae a ee 





where the denominators on the right are factors of @(X). The class 
found that, although this is not a sensible method to use, it works when 
the Y —r, are distinct as follows: 


f(x) = F(x) d(x) + R(x). 


For x=r;, (x) =0, therefore f(r;)=R(r,). The R(r,) is the number 
obtained on the left in the usual procedure of substituting 7; for x to 
obtain Ay, A2,-*:A,, Obviously, when ¢ contains repeated factors the 


method fails as we must substitute numbers other than r;. 


American Internationa! College. H. E. Bowler. 


Professor Albert E. Meder, Jr. tells us that the following song was 
sung to the tune of ‘The Lorelei” at a party of the Mathematics Club 
of the New Jersey College for Women, and that it is a translat‘on of a 
German poem by K. Lasswitz, 1880, which appears on pages 27 and 28 
of Leitzmann’s Lustiges und Merkwurdiges von Zahlen und Formen, 
Breslau, 1928. 
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SAD BALLAD OF THE JEALOUS CONES 


I. There once were two cones closely plighted 

In brotherly love and in bliss 

Upon the same axis united 
Their vertices joined in a kiss 

And then a smooth plane came a flying 
(Its smoothness was chief of its charms) 

And quickly, without even trying 
Seized one of the cones in its arms. 


II. An element chose at its pleasure— 
Then parallel took up its place 
While joy settled down without measure 
Upon the parabola’s face. 
The vertical cone found it tragic 
For never the gracious smooth plane 
Could cut it—not even by magic— 
It suffered in agonized pain. 


III. It urged and cajoled without ceasing 
And encouraged the plane soon to come 
With coaxing and begging and teasing 
It promised all comforts of home. 
The plane, which enjoyed being flattered 
Listened and heeded, entranced; 
A hyperbola now only mattered 
For beauty would thus be enhanced. 


IV. Serenely, with great satisfaction, 

The plane intersected each cone, 

But these were annoyed by its action, 
Each jealously deemed it its own 

Each jabbed with its point at its brother 
And tore at the unhappy plane 

On their axis they turned round each other 
And struggled with might and with main. 











Vi. 
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The plane cried, ‘‘Alas, cease your railings 
You're turning me round and round 
And though I despise finite failings 
I'll become an ellipse I'll be bound.”’ 


The cones, tho’ their struggle was needless 


Would cease not, not list to its call; 


And they pulled at each other, still heedless 


Of how the poor conic grew small. 


As the angry cones still madly battled 


A terrible shriek rent the air; 


The plane thru the vertex had rattled— 


The section had gone now for fair. 


Two surfaces battered and rumpled 


Forever and ever apart 


A crushed little point, badly crumpled 


The tragedy’s end breaks my heart!! 
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